NASA Technical Paper 1838 


NASA 
TP 
1838 
c . 1 



LOAN COPYo ^^S 
AF\VL TEGHa = 
KIRTLAND&^= 

t-* 

o 


A Microprocessor-Based Table 
Lookup Approach for Magnetic 
Bearing Linearization 


Nelson J. Groom and James B. Miller 


MAY 1981 




TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



□□bBiai 


NASA Technical Paper 1838 


A Microprocessor-Based Table 
Lookup Approach for Magnetic 
Bearing Linearization 


Nelson J. Groom and James B. Miller 
Langley Research Center 
Hampton, Virginia 


IWNSA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Branch 


1981 




SUMMARY 


An approach for producing a linear transfer characteristic between force 
command and force output of a magnetic bearing actuator without flux biasing 
is presented. The approach is microprocessor based and uses a table lookup to 
generate drive signals for the magnetic bearing power driver. An experimental 
test setup used to demonstrate the feasibility of the approach is described, 
and test results are presented. The test setup contains bearing elements simi- 
lar to those used in a laboratory model annular momentum control device (AMCD) . 


. INTRODUCTION 

This paper describes an approach for producing a linear transfer charac- 
teristic between the force command and force output of a magnetic bearing 
actuator. The approach, which is microprocessor based and uses a table lookup 
to generate drive signals for the magnetic bearing actuator power driver, was 
investigated for application to a laboratory model annular momentum control 
device (AMCD). The laboratory model (described in ref. 1) was built to inves- 
tigate potential problem areas in implementing the AMCD concept and is being 
used as part of an AMCD hardware technology development program. The basic 
AMCD concept is that of a rotating annular rim, suspended by a minimum of 
three magnetic bearing suspension stations and driven by a noncontacting elec- 
tromagnetic spin motor. A detailed discussion of the rationale for the AMCD 
configuration and of scxne of its potential applications is presented in 
reference 2. 

As described in reference 1 , the magnetic-bearing linearization technique 
used in the original laboratory model AMCD magnetic suspension was to differen- 
tially control sets of magnetic bearing elements about a permanent-magnet bias 
flux. Preliminary tests indicated that this approach (permanent-magnet flux 
biasing) presented a problem from a control system standpoint (ref. 3) . For 
a given equivalent permanent-magnet stiffness, a minimum bearing servo band- 
width is required for stability. The existence of structural modes in the area 
of the laboratory model bearing servo crossover restricted the amount of bear- 
ing servo damping that could be achieved. 

Because of the limitations of permanent-magnet flux biasing encountered 
with the laboratory model AMCD, a decision was made to explore an alternate 
approach to the design of the magnetic suspension system (ref. 4) . As a result, 
a new magnetic suspension system for the laboratory model has been designed, 
fabricated, and tested (ref. 5) . The new system uses a zero bias flux approach 
for the magnetic bearing actuators. Analog multiplier and square root modules 
produce a direct solution to the ideal magnetic actuator force equation to 
provide a linear force-current characteristic. (For further discussion of magr- 
netic bearing actuator control approaches, see ref. 6.) The accuracy of this 
approach is limited by the accuracy with which the ideal force equation approxi- 
mates the actual characteristics of the actuator and by the accuracy of the 
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analog components. This paper presents a zero bias flux linearization approach 
which is digital and which uses a lookup table constructed from measured actua- 
tor characteristics. An experimental test setup that was used to develop this 
approach is described, and test results are presented. 

SYMBOLS 

Dimensional quantities are presented in both SI Units and U.S. Custcxnary 
Units. Measurements were made in U.S. Customary Units. 

Fg force produced by bottom electrcxnagnet 

force command for magnetic bearing actuator 

F^(m) value of force command associated with mth line segment 

Fqi force produced by top electromagnet 

fg microcanputer system sample rate 

AFq = F^(m+1) - F^im) for all m, where m < N 

6Fc = Fc - Fc(m), where F^im) S F^ < F^Cm+l) 

G displacement of suspended element with respect to centered position 

in magnetic bearing actuator gaps 

Gg gap of bottom electromagnet 

Gq magnetic bearing actuator gap with suspended element centered 

Gip gap of top element 

Iq current in bottom electromagnet 

current command for magnetic bearing actuator 
Ic(m) stored value of current command associated with mth line segment 

current in top electromagnet 
K electromagnet constant 

N number of line segments 

SLOPE(m) slope of mth line segment from (F^ (m) ,1^ (m) ) to (P^ (m+1 ) ,1^ (m+1 ) ) 
Abbreviations: 

A/D analog-to-digital converter 
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AMCD cinnular mcmentum control device 

D/A digital-to-analog converter 

dc direct current 

emf electromotive force 


APPROACH 

Magnetic Bearing Control Approach 

The magnetic bearing control approach is one which uses zero bias flux. 
Figure 1 , a schematic representation of a magnetic bearing element pair, is 



i 

Go 
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Figure 1.- Magnetic bearing actuator. 


presented in order to describe this approach. Included in the figure are top 
and bottOTi electromagnets, with currents lip and Ig, respectively; a portion 
of the suspended element which is centered between the electromagnets; and a 
position sensor which measures the displacement G of the suspended element 
with respect to the centered position Gq. In the zero bias flux control 
approach, one electromagnet at a time is controlled; the top electromagnet is 
controlled for an upward force and the bottcxn electromagnet is controlled for 
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a downward force (since each electromagnet can physically produce only a uni- 
directional force) . Figure 1 indicates that if up is taken as the positive 
direction, the electromagnet gaps are 

Gt = Go - G (1 ) 

and 

Gg = Go + G (2) 

Assuming negligible fringing and ignoring nonlinear core effects, the force pro 
duced by each electromagnet is given by (ref. 4) 



The composite force-current characteristic of a zero bias flux actuator with 
the suspended element centered is shown in figure 2. 



Figure 2.- Composite force-current characteristic of a zero 
bias flux magnetic actuator- 



Linearization Approach 

The linearization approach is microprocessor based and uses a lookup table. 
The most straightforward approach for producing a linear transfer characteristic 
between force command and force output would be to solve equations (3) and (4) 
for the electromagnet current required to produce the desired force. That is. 


It = ^ 


Vcl\V2 


and 


Ib = Gg 


|Fcl\l/2 


(Fc > 0) (5) 


(Fc < 0) (6) 


where is the commanded force. Implementing the solutions of equations (5) 

and (6) digitally could produce a more computationally accurate result than the 
analog solution but would not account for deviation of the actual bearings from 
the ideal model. This solution would also require that considerable computing 
power be dedicated to a relatively minor portion of the total control system. 

The table lookup approach, which employs a one-to-one correspondence 
between force command and output current command permits the output/ 

input function to conform, within the limits of quantization error, to the 
actual bearing characteristics. However, this form of table requires con- 
siderable memory space. The minimum memory requirement is obtained by select- 
ing the minimum set of straight line segments which approximate the output/ 
input relationship within the desired tolerance. The major disadvantage of 
this choice of line segments is that considerable time may be spent searching 
for the appropriate line segment. 

Since in this application computational time is more critical than memory 
size, an approximation was selected which requires more memory but eliminates 
the search time. The minimum set of equally spaced line segments N are selec- 
ted so that N = 2*^, and the N + 1 points defining the line segments are sepa- 
rated by equal force command steps Af^. If the range of F^ is scaled to a 
B bit binary word (B > n) , the n most significant bits of F^ uniquely 
identify the N line segments. In actual computation these bits identify the 
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lookup table data associated with the mth force command data point F^^(m) such 
that Fc(m) ^ < F^ (m+1 ) . The B - n least significant bits of F^ repre- 

sent the difference 6 f^ between F^ and F^{m) . (See fig. 3.) 



Figure 3.- Table lookup approach. 


The data contained in the lookup table consist of the slope SLOPE (m) of 
each line segment m between (F^-, (m) , (m) ) and (F^;, (m+1 ) , (m+1 ) ) and of the 

current command Ic(rn) corresponding to F^(m). The location of the appropri- 
ate data is obtained directly from the force command word by masking and shift- 
ing. Calculation of output current command at G = 0 requires a single 

multiplication of 6 f^ by SLOPE(m) and addition of this result to I^Cm). 
Since actual variation of electromagnet current with gap was observed to be 
linear, as indicated by equations (5) and (6), calculation of for other 

values of G is accomplished by multiplying Iq by the appropriate bearing 
gap. For a more detailed description of the table lookup algorithm, see 
appendix A. 
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HARDWARE TESTS 


Description of Test System 

The test system consisted of a magnetic bearing test fixture connected to 
a microcomputer system as shown in figure 4. This system was used to obtain 
the data required to develop the lookup table and to obtain data on the perfor- 
mance of the proposed approach. A description of the current driver shown in 
figure 4 is given in appendix B. 
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Figure 4*- Magnetic bearing test system. 


Magnetic bearing test fixture .- The magnetic bearing test fixture shown 
in figure 5 consists of a magnetic bearing element pair; an equivalent 
"suspended" element, which is connected through a pair of load cells to the base 
of the fixture; and a pair of position sensors. The suspended element can be set 
to any desired vertical position in the magnetic bearing gap using the adjusting 
screws mounted on the load cells. The position sensors are used to measure the 
bearing gaps. The magnetic bearing elements have the same dimensions as the 
original magnetic bearing elements delivered with the laboratory model AMCD 
which is described in reference 1 . Two main differences exist between the test 
fixture bearing elements and the original laboratory model elements: (1) the 

core material of the test fixture bearing elements is SAE 1010 soft steel (as 
opposed to a lower loss silicon core iron used in the original elements) , and 
(2) the test fixture bearing elements contain no permanent-magnet material. 
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Figure 5.- Magnetic bearing test fixture. 


The load cells, shown in figure 5, are strain-gage bridge instrumented 
bending beams. The output of the bridge is a voltage which is directly pro- 
portional to the load applied to the beam. The cells were connected as shown 
in figure 6. They have a load range of ±44.48 N (±10 lb) and a nominal scale 
factor of ±0.045 mV/N-V (±0.2 mV/lb-V) . Scale factor and offset differ frcan 
cell to cell and vary with changes in test fixture configuration, power supply 
voltage, and temperature; therefore, software was designed to provide for 
periodic system calibration. 

Calibration was accomplished by applying a sequence of known loads to 
each cell and performing a first-order least-squares fit to the resulting data. 
Typical raw calibration data for one cell are given in figure 7. 
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Microcomputer system * - The same microcomputer system was used for both 
characterization and performance tests. The major system components include 
(1) microcomputer and memory, (2) analog interface unit, (3) cartridge data 
recorder, and (4) portable data terminal (fig. 4). 

The microcomputer is an Intel^ 8080 microprocessor-based computer which 
includes serial and parallel ports, programmable hardware timers, and 2000 bytes 
of random access memory (RAM) . This processor was selected for the initial 
phase of development because it is well supported in hardware and software. 
Available software includes FORTRAN as well as assembly language and the manu- 
facturer's high-level language. 

Most of the test programs were written in FORTRAN with assembly language 
hardware drivers. However, the table lookup routine was written in assembly 
language because of its time-critical nature. A listing of the table lookup 
routine is presented in appendix A. Memory expansion boards containing 24 000 
bytes of RAM memory were added to permit run time storage of large FORTRAN 
programs and substantial test data. 

The analog interface unit provides up to 32 multiplexed A/D input channels 
and two D/A output channels. The input channels provide the processor with 
position and force data from the test fixture position sensors and load cells 
and provide external command signals during those tests in which an analog force 
command is used. The two output channels are used to supply the current command 
to the bearing current driver. The analog interface unit was programmed to per- 
form conversions in two modes. During characterization tests, conversions were 
initiated under program control, and a flag on the interface was set at end of 
conversion. During those portions of the performance test when the table lookup 
algorithm was running in a real-time environment, the interface was programmed 
to start a conversion sequence on the rising edge of a hardware timer and to 
interrupt the processor at the end of conversion. This technique provides a 
stable sampling rate and does not require software timing loops. 

The data recorder was used for temporary storage and for transport of 
test programs and test results. This recorder permitted flexible use of the 
microcomputer system at a location which was remote fron the microprocessor 
development system used for software generation. 

The data terminal provided run time parameter selection, control, and 
monitoring of the system test programs. 


Description of Tests 

Characterization tests .- Characterization tests were required to establish 
the values of lookup table data. These data were collected by the microproces- 
sor system and transferred to a programmable desk-top calculator for analysis 


^Use of names of manufacturers in this report does not constitute an 
official endorsement of such manufacturers, either expressed or implied, by 
the National Aeronautics and Space Administration. 
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and reduction. A FORTRAN test program applied a 2048~point sequence of equally 
spaced current drive signals to the magnetic bearing test fixture and measured 
the actual force registered by the test fixture load cells at a manually selec- 
ted rim position. This process was repeated for a series of rim positions. 

These positions were varied frcmi -0.127 cm (-0.050 in.) to 0,1 27 cm (0.050 in.) 
in increments of 0.0127 cm (0.005 in.). The resulting force vs current data 
were converted by linear interpolation to a 2048-point current vs force table 
with equal force steps. These data were systematically reduced to form an 
N-point minimum size lookup table (N = 2*^) capable of approximating the measured 
data to within ±1 percent. This tolerance is within the limits of the existing 
analog solution to the ideal force equation. 

Performance tests .- Two tests were performed to evaluate the assembly 
language interpolation routine and the table lookup technique: 

1 . Linearity tests (static) - A sequence of force commands were applied 
to the test fixture through the linearization algorithm by a FORTRAN test pro- 
gram. The sequence of force commands was applied in the same direction and over 
the same range as in the characterization tests. The force produced on the rim 
was measured by the load cells and ccxnpared with the input command. 

2. Frequency response tests (dynamic) - A sinusoidal force-command signal 
was applied to the test fixture through the linearization algorithm, and the 
output current response was observed. The response measurements were performed 
by a frequency response analyzer. The force-command signal was given a dc 
offset so that measurements could be performed independently on upper and lower 
bearing elements. For this portion of the test, the characterization data were 
replaced by a linear table to permit direct comparison of input force command 
with current output. This replacement was necessitated by a mechanical resonance 
of the test fixture which prevented use of load cells for frequency response 
measurements. The frequency was measured over a bandwidth slightly larger than 
half the sampling frequency of the algorithm. 


TEST RESULTS AND DISCUSSION 
Characterization Tests 

The force, current, and gap relationships obtained from the data collected 
during characterization tests are similar to those obtained from the ideal 
electromagnet force equations (eqs. (3) and (4)). However, measured data 
deviated from ideal relationships near the zero force points, and differences 
between constants for upper and lower bearing elements were observed. Figure 8 
is a plot of the actual force data resulting from the application of a sequence 
of equally spaced current drive signals. Each curve represents 2048 data points 
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BERRING FDRCECN) 

Figure' 8.- Force-current data taken from magnetic bearing 

test fixture. 


taken at a given gap setting. The force vs current data with equal current 
steps were converted by linear interpolation to current vs force data with equal 
force steps. Analysis of the resulting data indicated that the relationship 
between current and gap for a given force was sufficiently linear to permit gap 
compensation to be accomplished by multiplying the output current by the gap (as 
in eqs. (5) and (6)). Systematic reduction of the data for G = 0 resulted 
in a 128-point table that was capable of reproducing the original data to within 
0.5 percent, as illustrated in figure 9. System memory requirements for this 
size table are quite reasonable. 


Performance Tests 

Linearity .- The main results of the table lookup algorithm linearity tests 
are shown in figures 10 and 11. Figure 10 shows the relation between force 
command input and the measured force output of the system. This particular plot 
is for G = 0, but the same result was obtained for other values of G. The 
actual percentage force error can be seen in figure 11. These results, although 
still within acceptable limits (less than ±1 percent error) , do not agree com- 
pletely with the errors predicted during data reduction (shown in fig. 9). 
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Figure 11 Percentage force error of table lookup algorithm* 


In particular, relatively large errors occur at the maximum positive and negative 
excursions of the force command. Since no attempt was made to minimize the 
effects of hysteresis in the test fixture, considerable residual magnetic flux 
can exist in the electromagnet cores and in the equivalent "suspended" element. 
Variations in the magnetic history of the test fixture are the probable cause 
of these errors. One possible way to reduce the effects of hysteresis on 
accuracy would be to use flux feedback in the power driver loop instead of cur- 
rent feedback. Another more obvious approach would be to use low-hysteresis 
material in both the electromagnet core and the suspended element magnetic cir- 
cuit material. 

Frequency response .- The table lookup algorithm required approximately a 
10-percent greater execution time than the 1.85-ms worst-case value predicted 
by evaluating the microprocessor instruction execution times. This loss of 
time occurred during the test because the table lookup algorithm instructions 
and data were stored in RAM memory on a memory expansion board. The increased 
execution time limited the system sample rate fg to a maximum of approxi- 
mately 490 Hz rather than the 540 Hz which was predicted. 

Results of the frequency response test are summarized in figure 12. The 
bandwidth exceeds 100 Hz, which was considered adequate for this portion of the 
bearing system. The magnitude is flat from dc to approximately fs/^ arid then 
rolls off gradually to approximately fs/^r which is the theoretical frequency 
limit for a sampled data system. Since this algorithm retains no history of 
the signals and has an almost fixed execution time, the phase response varies 
linearly with frequency to approximately fs/2* 
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CONCLUDING REMARKS 


A microprocessor based table lookup approach for magnetic bearing lineari- 
zation without flux biasing has been presented, and an experimental test setup 
used to demonstrate the feasibility of the concept has been described. Results 
obtained with the experimental test setup generally showed very close agreement 
with theoretical predictions. Using a 128-point table, the table lookup algo- 
rithm produced a linear transfer characteristic between force command and force 
output of the test fixture magnetic bearing actuator to within ±1 -percent error. 
The frequency response of the algorithm was greater than 100 Hz, which should 
be adequate for this portion of the bearing system. 

This approach when used as an inner loop for the actuators could form the 
basis for an all-digital magnetic suspension control system. One application 
would be the laboratory model annular momentum control device (AMCD). An all- 
digital system would allow control-system parameter changes to be made in soft- 
ware without requiring the circuit component changes and circuit rewiring which 
are necessary with existing analog systems. Also, advanced controller design 
approaches could be more easily implemented. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 19, 1981 
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APPENDIX A 


TABLE LOOKUP ALGORITHM 
Algorithm Description 

This appendix presents a listing, in assembly language, of the table look- 
up algorithm used to drive the magnetic bearing test fixture. The algorithm 
performs a nonlinear transformation of an input force command and a mea- 

sured rim displacement G into current commands I^, to either the top or 
bottom bearing element. The transformation approximates the solutions of equa- 
tions (5) and ( 6 ) in the main text. The lookup table contains 128 four-byte 
data groups. Each of the groups consists of a 16-bit value of SLOPE (m) and 
a 16-bit value of current I^Cni) obtained from measurement of the actual bear- 
ing force-current characteristic with the rim centered (i.e., Gqi = Gq - Gq) . 

The routine is interrupt driven and is called upon completion of the force 
command A/D conversion which is initiated by a hardware real-time clock. After 
reading the force command, the processor sets up the analog interface and initi- 
ates conversion of the rim displacement. Because of the configuration of the 
analog interface, the force command and rim displacement G are 1 2-bit 

two's complement numbers which are sign extended to 1 6 bits. 

The least significant 5 bits of F^ are removed and stored, since they 
contain the value The remaining 7 bits identify the 128 data groups. 

These bits are appropriately shifted, converted to offset binary, and added to 
the lookup- table base address to produce a memory pointer. The stored slope 
value is loaded using this memory pointer, and the pointer is incremented. The 
slope value is then multiplied by 6 f^ (an 8 -bit by 1 6 -bit multiplication rou- 
tine is used to conserve time) . This result is added to the stored value of 
l^(m) r which is now addressed by the memory pointer. The result, lc(^c'^o)' 
is the current command required to produce a desired force by either bearing 
element when the rim is centered (Gij» = Gb ” ^o^ • Since the required current 
is directly proportional to the bearing gap, the current command for any rim 
displacement is given by the following relationships: 


It - 



dc > 0) 


Ib - 



dc < 0) 


where Ic = f(Fc,Go)- After substitution of the table lookup algorithm for 
Iq and the expressions for Gip and G 3 from equations ( 1 ) and ( 2 ) in the 
main text, these equations become 
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It = ^lc(m) + SLOPE (m) 


Go - G' 




Ib = flc(i^) + SLOPE (m) 


(Go + G) 


Since Gq is a constant, these equations can be rewritten as 


It = 


Ic(m)\ /slope (m)\ 


I Go 


6f, 


(Go - G) 


- 


Ic(ro)\ /slope (m)\ 

+ 6f, 


(Go + G) 


The lookup table slope and current entries are predivided by G^ to eliminate 
a real-time multiplication by 1 /Gq. 

Because of the scaling of inputs and outputs (see table Al ) and the need 
to perform all calculations on 16-bit or less integer numbers, the stored look- 
up table values are scaled as shown in table A2. 


TABLE Al SCALING OF INPUTS AND OUTPUTS 


Variable 

Range 

Analog 
scale factor 

Analog 

range 

Digital 
scale factor 

Digital 

range 

Total 

scale factor 

Fc 

±11.12 N 
±2.5 lb 

0.9 V/N 
4 V/lb 

+10 V 

204.8 

±2047 

1 84.2 
81 9.2 

Ic 

±5 A 

1 V/A 

±5 V 

409.6 

±2047 

409.6 

G 

±0 .1 27 cm 
±0.05 in. 

378.7 V/can 
200 V/in. 

+ 1 0 V 

t>409.6 

±4096 

32.2 X 1q3 

5 X 2'' 4 


^These values include the sum of both position sensors. 

^his value includes a software multiplication by 2 (line 95) . 
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TABLE A2.- SCALING OF TABLE VALUE 


Variable 

Units 

Scaling factors^' 

SLOPE (m) 

A/K 

(0.5) (28) (2'' 6) 
(Go) (32.2) (1 03) 


A/lb 

(0.5) (28) (2'' 8) 
(Go) (5) (2^^!) 

Ic(m) 

A 

(2'' 6) 

(Go) (32.2) (2^4) 


A 

6) 

(Go) (5) (2^ 4) 


^The (0.5) factor accounts for the difference 
between the total scale factors for and 

^The (2®) and (2^^) factors are employed to maxi- 
mize the number of significant bits while maintaining 
intermediate and final results within a 16-bit integer 
format. These factors are removed by implied shifts 
during the two multiplications. 

^Factors (Gq) ( 32 . 2) (1 0 ^) and (Gq) {5) {2^ ate for 
Gq in SI (cm) and U.S. Customary Units, respectively. 




APPENDIX A 
Program Listing 

RSH88 :F1;L00KUP.RSH PRIHT(:F1;L00KUP.LST) 0BJECT( : FI: LOOKUP. 0BJ)PflGEHIDTH(88> 
EJECT 


ISIS-Il 8080/8085 MflCRO fiSSEMBLER, V3. 0 MODULE PAGE 1 


LX OBJ LINE SOURCE STATEMENT 




1 i 

3 i 


2 

^ LOOKUP. asM 



4 i 

5 

EXTRN CURPOS. CURNEG, ADCIN. TABLE. CONV. MUKADR. STAT 



US 

6 



^ SAVE MR IN PROGRAM REGISTERS 

0000 F5 


8 INTRU; 

PUSH PSH 


;SAVE A 4 PSH 

0001 C5 


9 

PUSH B 


. SAVE BC 

0602 D5 


10 

PUSH D 


.SAVE DE 

0002 E5 


11 

PUSH H 


iSAVE HL 



12 


13: 

; LOAD FORCE COMMAND AND 





14 

^ SETUP FOR GAP MEASUREMENT 

qocla oauuMi 
Cl WWW 

E 

15 

LHLD RDCIN 

;L0AD FORCE(N) 

0007 3E06 


16 

MVI A.6H 


>L0AD GAP CHANNEL 

0009 320000 

E 

17 

STA MUXADR 

; SELECT GAP CHANNEL 

000C 320000 

E 

18 

STA 


; START GAP CONVERSION 



19 


20 

> GENERATE DELTA FORCE COMMAND 





21 

> AND TABLE ADDRESS OFFSET 

000F 0600 


22 

MVI B.00H 


. CLEAR B 

0011 7C 


23 

MOV R.H 


iHIGH FORCE DATA > A 

0012 E60F 


24 

ANI 0FH 


iMRSK XT HIGH 4 BITS 

0014 67 


25 

MOV H.A 


iSTORE DBll THRU DBS > H 

0015 7D 


26 

MOV fl.L 


;L0H FORCE DATA > A 

0016 E61F 


27 

ANI IFH 


;MASK OUT HIGH 3 BITS 

0018 4F 


28 

MOV C.A 


; DELTA FORCE > C 

0019 fiD 


29 

XRA L 


>M0VE LOW DATA > A (LOH 5 BITS M 



ASKED) 




001fl 84 


30 

ADD H 


.•COMBINE DBll THRU DBS WITH DB7 



THRU DBS 



001B 07 


31 

RLC 


} 

001C 07 


32 

RLC 


^ REARRANGE BVTE 

0010 07 


33 

RLC 


; TO FORM 

001E 07 


34 

0 

35 

RLC 


; DB10. DB9. DBS. DB7. DB6. DBS. 0. 

001F 17 


RAL 


i 

0020 6F 


36 

MOV LA 


; CONVERT TO 

0621 3F 


37 

CMC 


; OFFSET 

0022 17 


38 

RAL 


; BINARY 

0023 E601 


39 

ANI 01H 


; FORM 

0025 67 


40 

MOV H.A 


; 



41 


42 

^ COMPUTE ADDRESS OF SLOPE <N> 

0026 110000 

E 

43 

LXI D. TABLE 

;L0A0 TABLE START ADDRESS 

0029 23 


44 

INX H 


; ADJUST nABLE" ADDRESS VALUE 

002fl 23 


45 

INX H 


} 

002B 19 


46 

DAD D 


>ADD TABLE START TO OFFSET 



47 






48 


49 

> LORD SLOPE CN> 
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ISIS-II 8888/8885 NflCfiO fiSSEWLER, V3, 8 NOOULE PflGE Z 


LOC OBJ 


LINE 

SOURCE ST8TENENT 


002C 5E 


50 

H0VE.H 

LOflO LOW BYTE SLOPE(N) 

0020 23 


51 

INX H 

INCREMENT POINTER FOR HIGH BYTE 

002E 56 


52 

MOV 0,M 

LORD HIGH BYTE SLOPE(N) 

002F 23 


53 

INX H 

INCREMENT POINTER FOR CURRENT fl 



DORESS 





54 

55 

? MULT I PL V SLOPE CN> BV DELTA FORCE 




56 

^ <SPECIAL SBIT X 16BIT MULT. > 

0038 E5 


57 

PUSH H 

SAVE ADDRESS 

0031 EB 


58 

XCHG 

SHflP MULTIPLICflN(SLOPE) TO HL 

0032 220000 


59 

SHLD TEHP 

SAVE MULTIPLICflN IN TEMP 

0035 21CE00 


68 

LXl H.BNUH 

LOAD CVaE COUNTER ADDRESS 

0038 3609 


61 

MVI H.09H 

LOAD CVaE COUNTER 

0038 110000 


62 

LXI D.0H 

CLEAR TEMPORARY RESULT 

0030 79 


63 LOOP: 

MOV fl;C 


003E IF 


64 

R8R 

ROTATE MULTIPLIER<DELTAF> 

003F 4F 


65 

MOV C.R 


0040 35 


66 

OCR M 

DECREMENT CYCLE COUNTER 

0041 CR5800 


67 

JZ FINI 

TEST FOR MULTIPLY COMPLETE 

0044 O24F00 


68 

JHC SKIP 

JUMP IF MULTIPLIER BIT 6 

0047 2BO000 


69 

LHLO TEMP 

GET MULTIPLICflN 

004fl 19 


70 

080 0 

AND ADD 

004B EB 


71 

XCHG 

SAVE PARTIAL PRODUCT 

004C 21CE00 


72 

LXI H,BNUM 

RELOAD BNUM ADDRESS 

004F Tfl 


73 SKIP: 

MOV 8;D 


0050 IF 


74 

RRR 

ROTATE 

0051 57 


75 

MOV 0,8 


0052 7B 


76 

MOV fl,E 

TEMPORARY 

0053 IF 


77 

RRR 


0054 5F 


78 

MOV E,fl 

RESULT 

0055 C33D00 


79 

JMP LOOP 

LOOP 

0058 El 


80 FINI: 

81 

POP H 

RESTORE ADDRESS 



82 

; COMPUTE CURRENT COMMAND FROM 




S3 

; BASEICN> AND DELTA I 

0059 4E 


84 

MOV C,M 

LOAD LOW BYTE BflSEKN) 

0058 23 


85 

INX H 


005B 46 


86 

MOV B>M 

LOAD HIGH BYTE BASEKN) 

005C EB 


87 

XCHG 

SWAP DELTAI TO HL 

0050 09 


88 

DflD B 

;ADD BASEKN) TO DELTAI 

005E 44 


89 

MOV B,H 

iMOVE CURRENT COMMAND 

005F 40 


90 

91 

MOV C/L 

i TO BC 



82 

; CORRECT FOR BEARING GAP VARIATION 

ww jfiCTJw 

E 

93 

LOB ST8TUS 

laEAR RTI1200 EOC FLAG 

0863 280008 

E 

94 

LHLO fiOCIN 

iLOAD GAP 

0066 29 


95 

DflD H 

.DOUBLE GAP 

0067 3E05 


% 

MVI asH 

;LOAD FORCE COMMAND CHANNEL 

0069 320800 

E 

97 

57fl WXflDR 

: SELECT FORCE COMMAND CHANNEL 

006C 111621 


98 

LXI D,GflP 

iLOAD NOMINAL GAP VALUE 

006F B0 


99 

IS 

ORB B 

1 TEST FOR SIGN OF CURRENT .COMMAN 

0070 F5 


U 

100 

PUSH PSH 

; SAVE SIGN FLAG 

0871 F27F00 


101 

JP UPPER ^JUMP TO UPPER BEARING COIL 



102 




21 



APPENDIX A 


ISIS-II 8080/^ rifiCRO ASSEMBLER. VI 0 

MODULE PAGE 3 


Loc oej 

LINE 

SOURCE STATEMENT 





±03 

j COMPUTE LOWER BERRING 

COIL GRP 

0074 19 

IW LOWER: 

DAD D 

;ADD GAP(T) TO NOMINAL GAP 


0075 00 

165 

DCX B 



0076 79 

106 

MOV A.C 

;TAKE 


0077 2F 

10? 

CMA 

; ABSOLUTE 


0078 4F 

108 

MOV C,A 



0079 78 

109 

MOV A.B 

; VALLE 


007fl 2F 

U0 

CMA 

; OF 


0078 47 

111 

MOV B.A 

; CURRENT 


007C C38500 

112 

JMP FINAL 

;JUMP TO FINAL CURRENT COMMAND C 



OMPUTflTION 




112 





114 

±±5 

j COMPUTE UPPER BERRING 

COIL GRP 

007F 78 

116 UPPER: 

MOV A.E 

i SUBTRACT 


0080 95 

117 

SUB L 

i GAP<T) 


0081 6F 

118 

MOV LA 

; FROM 


0082 7ft 

119 

MOV A.D 

; NOMINAL 


0883 9C 

120 

SB8 H 

; GAP AT 


0084 67 

121 

MOV H.A 

; CENTER 



122 






±23 

; MULTI PLV GRP BV CURRENT COMMRND 

0085 220000 

124 FINAL; 

SHLD TEMP 

; STORE MULTIPLICAN IN TEMP 


0088 21CE00 

125 

LXI H.BNII1 

; STORE 


0088 3611 

126 

MVI M.11H; 

i BIT COUNT 


0080 110000 

127 

LXI D.00H 

; INITALIZE RESULT 


0090 78 

128 LOOPl: 

MOV A.B 

; ROTATE 


0091 IF 

129 

RAR 

f 


0092 47 

130 

MOV B/A 

; MULTIPLIER 


0093 79 

131 

MOV A.C 

/ 


0094 IF 

132 

RAR 

; RIGHT 


0095 4F 

133 

MOV C.A 



0096 35 

134 

DCR M 

; DECREMENT BIT COUNT 


0097 CftHE00 

135 

JZ FINIl 

jOONE? then output 


009ft O2ft500 

136 

JNC SKIPl 

; JUMP IF NO CARRY FROM ROTATE 


0090 2ftD000 

137 

LHLD TEMP 

; OTHERWISE 


0000 19 

138 

DAD D 

; ADD MULTIPLICAN 


00ftl EB 

139 

XCHG 

;SAVE RESULT 


00ft2 21CE00 

140 

LXI H.BNUM 

; RESTORE BIT COUNT POINTER 


00ft5 Tfl 

141 SKIPl: 

MOV A.D 

i ROTATE 


0086 IF 

142 

RAR 

; TEMP 


00fl7 57 

143 

MOV D.A 

; RESULT 


0008 78 

144 

MOV A.E 

; RIGHT 


0009 IF 

145 

RAR 

; 


0000 5F 

146 

MOV LA 

j 


0008 C39000 

147 

JMP LOOPl 

.REPEAT LOOP 



148 

±49 

; CHECK FOR OVERCURRENT 

VRLUE 



±50 

j RND OUTPUT RESULT TO CORRECT PORT 

000E EB 

151 FINIl: 

XCHG 

^SHAP CURRENT TO HL 


000F 7C 

152 

MOV A.H 

;LOAO CURRENT 


0080 FE08 

153 

CPI 08H 

;TEST CURRENT VALUE 


0082 Dft8800 

154 

JC OK 

jJUMP IF OK 


0085 21FF07 

155 

LXI H.07FFH 

;CURRENT=MRXIMUH CURRENT 


0088 FI 

156 OK: 

POP PSH 

i RESTORE CURRENT SIGN 
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APPENDIX A 


I5I5-II 8080/8885 HflCRO ASSEMBLER/ VI 0 


MODULE PAGE 


LOC 

OBJ 


LINE 

SOURCE STATEMENT 



0063 

F2C200 


157 

JP OUTPOS 

/JUMP IF POSITIVE 


006C 

220000 

E 

158 

SHLD CURNE6 

/OUTPUT RESULT TO 

LONER 




D/A 




00BF 

C3C500 


159 

JMP RETURN 

i RETURN 


00C2 

220000 

E 

168 OUTPOS: 

SHLD CURPOS 

; OUTPUT RESULT TO 

UPPER 




D/A 




00C5 

El 


161 RETURN: 

POP H 

/RESTORE HL 


00C6 

D1 


162 

POP D 

; RESTORE DE 


00C7 

Cl 


163 

POP B 

^RESTORE BC 


00C8 

3E20 


164 

MVI A/20H 

; LOAD END OF INTERUPT 

00CA 

D30A 


165 

OUT 0DAH 

; OUTPUT EOIC 


00CC 

FI 


166 

POP PSH 

/POP A A PSH 


00CD 

C9 


167 

RET 

; RETURN 





168 







169 







170 







171 




00CE 

QOM 

TOW 


172 BNUM: 

DN 0 



0000 

OwO 


173 TEMP: 

ON 0 



2116 



174 GAP 

EQU 8470 






175 

END 




PUBLIC SYMBOLS 


EXTERNAL SYMBOLS 



ADCIN E 0000 

CONV 

E 0000 

MUXAOR E 0000 

STATUS 

E 0000 

USER SYMBOLS 



ADCIN E 0000 

BNUM 

A 00CE 

CURPOS E 0000 

FINAL 

A 0085 

GAP A 2116 

INTRU 

Q fWV> 

n Uwv 

LONER A 0074 

MUXADR 

E 0000 

RETURN A 00C5 

SKIP 

A 004F 

TABLE E 0000 

TEMP 

A 00D0 

ASSEMBLY COMPLETE/ NO 

ERRORS 


CURNEQ E 0000 

CmOS E 0000 

TABLE 

E 0000 


CONV 

E 0000 

CURNEG E 0000 

FINI 

A 0058 

FINIl A 00AE 

LOOP 

A 003D 

LOOPl A 0090 

OK 

A 0068 

OUTPOS A 00C2 

SKIPl 

A 00A5 

STATUS E 0000 

UPPER 

A 007F 



APPENDIX B 


MAGNETIC BEARING CURRENT DRIVER 

A schematic diagram of the current driver used in the magnetic bearing 
test system is shown in figure B1 . The driver is capable of supplying up to 



Figure Bl Magnetic bearing current driver. 


5 A to each bearing coil and has a gain of 1 A/V. Drive is provided by the 
LH 0002 current amplifier at low current levels (<100 mA) and by the comple- 
mentary Darlington configuration at higher levels. The driver provides flat 
response over a bandwidth greater than 200 Hz (fig. B2) . The frequency 
response is limited by the ability of the power supply voltage to overcome the 
back emf of the bearing coil, and the effect is observed as an apparent slew 
rate limitation. The 200-Hz bandwidth is sufficient to evaluate the table look- 
up algorithm. 
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